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^ ' Abstract 

^ . We calculate the radiative corrections to the unpolarized and the four polarized 

Q I spectrum and rate functions in the leptonic decay of a polarized /i into a polarized 

' electron. The new feature of our calculation is that we keep the mass of the final 

^> . state electron finite which is mandatory if one wants to investigate the threshold 

^ I region of the decay. Analytical results are given for the energy spectrum and the 

^ ' polar angle distribution of the final state electron whose longitudinal and transverse 

polarization is calculated. We also provide analytical results on the integrated spec- 
trum functions. We analyze the nif, limit of our general results and investigate 
the quality of the mg — > approximation. In the mg — > case we discuss in some 
detail the role of the 0(a) anomalous helicity fiip contribution of the final electron 
which survives the me — limit. The results presented in this paper also apply to 
the leptonic decays of polarized r-leptons for which we provide numerical results. 



1 Introduction 



While calculating the radiative 0{as) QCD corrections to polarization effects in semilep- 
tonic decays of heavy quarks, where the full quark mass dependence was retained [1, 2], 
we came to realize that our results could also be gainfully employed in the corresponding 
0{a) QED corrections to the weak leptonic decays of the /i- and r-leptons^. In most of 
the previous radiative correction calculations the mass of the charged lepton daughter /' 
has been neglected except for anomalous contributions from the collinear region which 
survive the m'l —>■ limit [3, 4, 5] and the logarithmic terms ~ (Inm;) which are needed 
to regularize the collinear divergencies that appear in the loop and tree graph ("internal 
bremsstrahlung" ) contributions. These logarithmic terms partially cancel in the spectrum 
and completely cancel in the rate when the loop and tree graph contributions are added. 

From general considerations it follows that the unpolarized and three of the polarized 
spectrum functions contain only even powers of the mass ratio m^/m/. Considering the 
fact that {mjm^f = 2.34 ■ 10-^ {m^/m^f = 3.54 ■ 10^^ and {mjm^f = 8.27 • 10"* 
the zero mass approximation should be a good approximation for most of the energy 
spectrum of the daughter leptons except for the region close (or very close) to the soft 
endpoint of the spectrum (also referred to as the threshold region) where finite mass 
effects have to be retained. Contrary to this the transverse polarization of the daughter 
lepton is proportional to the linear mass ratio mj/ mi. Also, when integrating the spectrum 
functions, the linear mass ratio enters in all four polarized rate functions. Finite mass 
corrections may thus play an important role at least for (r —>■ /i)-decays where the linear 
mass ratio mn/mr = 5.95 ■ 10~^ is not very small. An improved analysis of r-decays is of 
quite some topical interest since large samples of r-leptons are currently being produced 
at the existing two B factories in Japan and in the USA, and are expected to be produced 
at future r-charm factories to be set up in Ithaca and Beijing. As the data becomes more 
precise, the predictions of the SM including also radiative correction effects will be tested 
at an ever-rising level of precision. It is then evident that the inclusion of final lepton 
mass effects must play an important role in particular in the threshold region. 

We determine the radiative corrections to the daughter's lepton energy spectrum and 
its longitudinal and transverse polarization keeping the dependence on the polarization 
of the parent lepton. This generalizes the calculation of [6] in which the zero mass ap- 
proximation was used. Our calculation extends the calculation of [7], which also includes 
finite mass effects, in that we include the longitudinal and transverse polarization of the 
daughter lepton. 

^The three leptonic ^- and r-decays are treated within the Standard Model, i.e. they are all governed 
by the same {V — A) coupling structure and coupling strength Gp- 
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The paper is structured as follows. In Sec. 2 we introduce our notation and write 
down the general structure of the spin-dependent rate. Sec. 3 contains our Born term 
results. ly-propagator effects are taken into account in Sec. 4. In Sec. 5 we present our 
analytical and numerical results on the 0{a) radiative corrections. In Sec. 6 we consider 
the m'l ^ limit of the mj 7^ results presented in Sec. 5 and discuss in some detail 
the origin of the anomalous helicity flip contribution. Sec. 7 contains our summary and 
conclusions. In two Appendices we collect some technical material on trilog functions and 
Fierz identities relevant to our calculation. 

2 General structure of spin-dependent rate 

To make life simple we shall in the following always refer to the specific case fi^ 
e~ + Ue + i^fi instead of referring to the generic case involving also leptonic r-decays when 
writing down analytical results. Of course, when discussing numerical results, the two 
leptonic r-decay channels r~ — /i~ + + u^- and r~ ^ e~ + + are also included. 

Since the subject of /^"-decays is extensively covered in the relevant textbooks (see 
e.g. Refs. [8, 9, 10]) and review articles (see e.g. Refs. [11, 12, 13]) we can afford to be 
very brief in describing the formalism. 

From helicity counting one knows that there are altogether five spin-dependent struc- 
ture functions and one spin-independent structure function describing the leptonic decay 
of a polarized muon into a polarized electron. We thus define a spin-dependent differen- 
tial rate in terms of six invariant structure functions Ai. In the rest system of the the 
decay distribution reads 



As usual X = 2Ee/m^ denotes the scaled energy of the electron where the energy of 
the electron is defined in the rest frame of the 9p is the polar angle between the 
polarization of the muon and the momentum direction of the electron in the muon rest 
frame. 

Eq. (1) will be evaluated in the rest system of the muon where = {m^\ 0, 0, 0) 
and Pe = {Ee;0,0,\pe\) = (m^/2)(x; 0, 0, x/3). The velocity of the electron is denoted 



dr 




dx dcos 6 




(1) 
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hy P = yl — Ay'^/x^ where y = me/m^. In the rest frame of the fx the polarization 
four-vectors of the fi~ and e~ are given by 

5^ = (0;C), (2) 

where the polarization three-vector (^^ of the /i~ and the quantization axis fie of the spin 
of the e~ in their respective rest frames read (see Fig. 1) 

4 = (sin^P, 0, cos^p) (4) 

and 



fie = (sin^cosx? sin 6^ sin X, cos^). (5) 

Eq. (4) holds for 100% polarized muons. For partially polarized muons with magnitude 
of polarization P the representation (4) has to be multiplied by P such that = PC^- The 
representation Eq. (5) needs a word of explanation. The vector fie denotes the orientation 
of the quantization axis of the electron's spin in the rest frame of the electron which can 
be freely chosen. The orientation of fie has no physical meaning per se. In particular fie is 
not the polarization vector Pe of the electron whose Cartesian components can, however, 
be obtained by projecting onto the a;-axis (6* = 7r/2, x = 0), the y-axis {9 = 7i/2,x = 7r/2) 
and the z-axis {6 = 0). 

One finally has [6, 14] 

- — f^—— = Px To{Gi + G2P cos 9p + Gs cos 9 + G^P cos 6p cos 9 
ax acos 9p 

+ G^P sin 9p sin 9 cos X + Ge-P sin ^p sin ^ sin x). (6) 



The relation between the invariant structure functions Ai and the frame-dependent 
spectrum functions Gi is given by 



Gi = A,, 
G2 = -^xpA2, 
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2y 



4y 2y 



G5 — —A^, 

Ge = \x(3A,. (7) 

Gi is the unpolarized spectrum function, G2 and G^ are single spin polarized spectrum 
functions referring to the spins of the /i~ and e~, resp., and G4, G5 and Gq describe spin- 
spin correlations between the spin vectors of the muon and electron^. Gq represents a so- 
called T-odd observable. This is evident when rewriting the angular factor multiplying Gq 
in (6) in triple-product form, i.e. sin 6'p sin ^ sin x = |pe|~^Pe ■ (Ca* ^ ''^e) ^- Gq is identically 
zero in the SM since, on the one hand, the weak coupling constant Gp real and, on the 
other hand, the loop contributions do not generate imaginary parts. Gq will therefore not 
be discussed any further in the following. 

It is quite instructive to rewrite Eq. (6) in a factorized form 

/5a;ro(Gi + G2Pcos^p)(l + A-ne) (8) 



dx dcos Op 

A comparison with Eq. (6) shows that the polarization vector P^. of the electron has the 
components 



py ■— pN 



e 



Gi + G2Pcosep' 

PGq sinOp 
G1 + G2P cos Op' 



I ^ Gs + G^P cos 9p 

Gi + G2Pcosep ^ ' 



^The spectrum function G5 describing the transverse polarization of the electron vanishes for vanishing 
electron mass and is therefore not included in the corresponding decay distribution in [6] . 

•^A general discussion of the electron polarization in muon decay can be found in [15], including a 
discussion of tests of the TCP theorem. 
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We have as usual denoted the (x, y, z) components of by {P;!~, P^, P^) where (±, A^, /) 
stand for the polarization components transverse to the momentum direction of the elec- 
tron (in the plane spanned by the momentum of the electron and the polarization vector 
of the muon), normal to this plane and longitudinal, respectively. Eq. (9) shows that 
the spectrum functions and G4 determine the longitudinal polarization of the electron 
while its tranverse polarization (in the plain spanned by the electron's momentum and 
the muon's polarization) is determined by Gs^. 

The limits on x = 2Ee/m^ are given by Xmin = 2?/ and Xmax = 1 + 1/^ where, as before, 
y = me/mfj,. Fq, finally, is the me = Born term rate given by Fq = G|^m^/1927r^. The 
differential rate for the charge conjugated decay /i^ — + z/e + i^ti is obtained from 
Eq. (6) by the substitution Gj — > Gj (i = 1, 4, 5, 6) and Gi —Gi {i = 2, 3) [15]. 

It is convenient to split the unpolarized and polarized rate functions into a Born term 
part and a radiatively corrected part according to 

Gi = Gf°™ + GS"^ z = 1,2, 3, 4, 5. (10) 

The respective results on the Born term contributions Gf°™ and the 0{a) corrections to 
the rate functions G^^^ are given in Sees. 3 and 5. The sum of the two contributions in 
Eq. (10) will generally be referred to as the next-to-leading order (NLO) result. 



3 Born term results 

We shall work with the charge retention form of the Lagrangian for the decay ^~ 
e~ + z/g + which reads^ 

C{x) = ^t,(x)7"(l - 75)vl>^(x)'I'.,Xx)7„(l - 75)*.. (x) + h.c. (11) 

When squaring the corresponding matrix element one obtains the tensor G"^ from the 
charged lepton side (G for charged) which has to be contracted with the neutrino-side 

*As noted before in the SM there is no transverse polarization normal to that plane. 

^The charge retention form of the Lagrangian is obtained from the usual Standard Model charged 
current-current form through the Fierz transformation (Bl) written down in Appendix B. The minus 
sign from the Fierz identity is cancelled from having to commute the Fermion fields an odd number of 
times in order to relate the two forms. Note that the Fierz identity is a four-dimensional identity. Since 
our calculations are done in four dimensions the Fierz identity can be safely applied. 

^When the Lagrangian for /^-decay is written in charge retention form the similarity of the decay 
e~ + Ve + to the decay 5 — > c + e~ + i>e becomes quite apparent through the substitutions 
b <-^- c <-> e^, i^fi and Pe f'e- 
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tensor N"^ {N for neutral). For the Born term contribution of the charge-side tensor one 
obtains 

Csin = \^^{il/>e + me){t + 75^)7^(1 ' l^){i>, + + 75^^)7^(1 " 75)}, (12) 

where the dependence on the polarization four-vectors of the /i~ and e~ has been retained. 

Since only even-numbered 7-matrix strings survive between the two (1 — 75)-factors 
in Eq. (12) one can compactly write the result of the trace evaluation as 

Cforn = nP^.P: + P^Pf - g'-^P.-pe + te'^^^' p,,,p,,s) , (13) 

where 

Pe=Pe-^e< (15) 

and where and Sg are the polarization four-vectors of the fi~ and e~, resp., defined in 
Eq. (2). 

The dependence on the momentum directions of the z/g- and z/^-neutrinos has been 
completely integrated out in the differential rate. Thus the neutrino-side of the interaction 
can only depend on the spatial piece of the second rank tensor build from the momentum 
transfer to the neutrinos (for the present purpose the neutrinos are treated as massless) 
which we denote by ''. Thus the relevant neutral-side tensor is given by 



Q 

The Born spectrum functions can then be extracted from 



= -g-^ + (16) 



Q^J^aP^Born ^ 2((p^-pe) + 2{p,.Q) {p^-Q)), (17) 

where the antisymmetric piece in C^p'^"' has dropped out after the symmetric contraction. 
At the Born term level one has Q = Pfi — Pe- 

Including the correct normalization the differential Born term rate is given by 



^We denote the momentum transfer to the neutrino pair by a capital Q in order to set it apart from 
the momentum transfer q to the (e~i>e)~pair used in Sec. 4. 
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Born 



dxdcosOp 

The spectrum functions defined in Eq. (6) can then easily be calculated from the Born 
term contributions (17) using the relations Eq. (7). They are given by 



^Born 


x(3- 


- 2x) - (4 - 3x)y^, 


^Born 


(3x{l 


-2a; + 32/2), 


^Born 


-px{3 


-2X + 2/2), 


QBorn 


-x{l - 


- 2x) - (4 + x)y'^, 


^Born 


-2y{l 


- X + y"^). 



(19) 



Note that the Born term spectrum functions G'f2™4 and G^°^"'/y are quadratically de- 
pendent on y. This is in agreement with the general arguments presented in [22]. For 
y^ = {rrie = 0) one has Gf °™ = -Gf and Gf °™ = -Gf °™ which reflects the fact 
that a mass zero left-chiral electron is purely lefthanded. The Born term results (19) 
reproduce the y = results of Ref. [6]. For y"^ ^ {jrie 7^ 0) our Born term results for 
QBorn Gf agree with those of Ref. [7] . Note that the spectrum function Gf is 
proportional to y and thus vanishes for vanishing electron mass. The overall chiral fac- 
tor y = mf./m^ in G5 originates from a lefthanded/righthanded interference contribution 
which is chirally suppressed. 

In Figs. 2a (yU ^ e), 3a (r fi) and 4a (r e) we show plots of the a;-dependence 
of the four Born term spectrum functions /3a:Gf2^^4. They rise and fall from zero at 
the soft end of the spectrum to (1 — y'^Y {i = 1,4) and —(1 — y'^Y {i = 2,3) at the 
hard end of the spectrum, respectively. The mi' = pattern Gf'"'" = — Gf°™ and 
QBorn _ _QBorn jg giig]-^^iy distorted by final lepton mass effects except for the point 
Xmax = 1 + where the above = relations are exact. For the spectrum function 
Gf ™ one finds Gf ^'^(x^^^ = 2y) = -2y{l - yf and Gf ™(x^a, = l + y'^) = Q. The 
chirally suppressed contribution of (3xG^°^^ is negative over the whole x-range. At the 
scale of the plots it is only visible for the case r — /i (Fig. 3a). 

In addition to the polarization observables we also define a forward-backward asym- 
metry for the case that the spin of the electron is not observed. It reads 

Fp — Fr 1 „ Go / 
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where Tp and Tb are the rates into the forward (cos^p > 0) and backward {cos 6p < 0) 
hemispheres. 

One can also define a forward-backward asymmetry of the longitudinal polarization 
proportional to G4/G1 according to 

where the indices +/— denote the helicities of the electron. This asymmetry will be diffi- 
cult to measure since it involves a spin-spin correlation measurement. We shall therefore 
not discuss this asymmetry any further in this paper. 

The longitudinal polarization and the forward-backward asymmetry take the values 
pi,Born _ _ippQg^^ Qj^^ App"^ = 0, respectively, at the lower limit Xmin = 2?/ where 
QBorn ^ _3G'iJo™ qBotu ^ qBovu ^ q rpj^j^ contrasted with the naive limits 

pi^Born _ g^^^ App^"' = |P when naively setting y = in the corresponding ratios. At 
the upper limit Xmax = (1 + 2/^), where Gf°™ = -Gf"™ = -Gf''™ = Gf°™ = (1 - y^f, 
the longitudinal polarization decreases to = — 1 irrespective of the value of P cos 6p, 

and the forward-backward asymmetry decreases to App"' = — |P as in the naive me = 
case. At the soft end of the spectrum one finds a substantial transverse polarization 
p±,Born _ — Ipsin^p which is independent of y contrary to naive expectations. Thus 

-» Born 

the total polarization of the electron at threshold is given by |Pe | = ^P irrespective 
of the value of cos6p. At the hard end of the spectrum one has p^^'^°™ = 0. 

In Figs. 2b (yU. — > e), 3b (r — > fi) and 4b (r — > e) we show the x-dependence of the Born 
term prediction for the longitudinal and transverse polarization of the daughter lepton. 
We set P = 1 and take three values cos6p = 1,0 and —1 for the longitudinal polarization 
and set cos6p = for the transverse polarization. The longitudinal polarization stays 
close to —1 over most of the (hard part) of the spectrum but deviates significantly from 
the naive value —1 in the threshold region with only a slight dependence on the value 
of cos 6'p. In order to highlight the deviations from the naive value p^';-^"''" = — 1 in 
the threshold region we have chosen logarithmic scales for the energy variable x. The 
transverse polarization is negative and stays very close to zero over most of the (hard 
part) of the spectrum and decreases to its limiting value —1/3 at threshold. In fact all 
the Born term curves can be seen to approach the limits discussed above at the soft and 
hard end of the spectrum. 

In Figs. 2c, 3c and 4c we show the corresponding curves for the forward-backward 
asymmetry A^'^^ for the three decay cases. Again we set P = 1. The forward-backward 
asymmetries rise from the limiting value A^g^" = at threshold and then fall to A^'^'^ = 
—1/2 at the hard end of the spectrum. 
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-0.019% 


-0.12% 


-2.82% 


- 8.48% 


-6.62 


10- 


-5% 


-7.39 ■ 10-4% 


G2 


-3.57-10-4% 


-0.67% 


-0.57% 


-11.23% 


-7.60 


10- 


-«% 


-0.039% 


G3 


-0.031% 


-1.48% 


-4.25% 


-24.49% 


-1.10 


10- 


-4% 


-0.084% 


64 


-0.019% 


-1.43% 


-2.82% 


-22.15% 


-6.62 


10- 




-0.083% 



Table 1: 

Percentage mass corrections to m// = (Born) and mi' — > (0(a)) rate functions for the three 
cases {fi ^ e), (r — > fj.) and (r — > e). 



Next we integrate the differential Born term rates over the full x spectrum. Let us 
define reduced Born rate functions Gf°^^ according to 

1+3/2 

Gf = jdx(3x Gf °™ i = 1, 2, 3, 4. (22) 
22/ 

y4)(l - 81/2 + ^4) _ i2y4in^ = i - 42/2 + 0(y4), 
yf{l + 5y + 15y2 + 3y3) = _ 1 + i^^3 + 0(^4)^ 

^)5(3 + + 5^2 + ^3) ^ _i + _ ig^3 + 0(^4)^ (23) 

y4)(l - 81/2 + ^4) _ 4^4 ij^^ ^ 1 _ 4^2 + 0(2/4), 
GBo™ ^ _^(i(i _ y2)(i _^ 10^2 + ^4) + 4^2(1 + ^2) 1^^) ^ + (3 + 4 l^^)^2 + ^(1/4)) 

The occurrence of odd powers of y in and (5^°™ can be traced to the lower 

boundary x = 2y oi the x-integration which is linear in y. In view of this it is quite 
remarkable that Gf°^", Gf°^"' and G^°^"'/y contain only even powers of y. Nevertheless, 
the leading mass correction to sets in only at 0{y^) and, for only at 0{y^). In 

this sense the Born term rates are well protected against finite electron mass effects, or, in 



One obtains 



^Born 


1 


:i 


^Born 




:i 


^Born 




:i 


/yBorn 
<-J4 — 


1 
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the case (t ^ /i) , reasonably well against muon mass effects. This is illustrated in Table 1 
where we list the values of the Born term percentage changes {Gi{m 7^ 0) — Gj(m = 0)) 
/Gi{m = 0) when going from m^/ = to mi' 7^ for all three cases. 

The average longitudinal polarization of the electron (Pj) and the average forward- 
backward asymmetry {Afb) is obtained by the replacement {Gi Gi) in (9) and (20). As 
in the rate expressions (23) the final state lepton mass effects on the Born term polarization 
(Pg) and forward-backward asymmetry {Aps) are quite small. The deviation from the 
mil = result (P//) = —1 is of O(10~^) in muon decay and O(10~^) in (r — > e). For 
(r ^ fi) the deviation from (P^) = —1 is of O(10~^). The dependence on the value of 
P cos 6p is very small. The average value of the forward-backward asymmetry (Afb) is 
quite close (0(10^^ — 10^'^)) to they = prediction (Afb) = —1/6. The average transverse 
polarization (P^) is generally quite small due to the overall chiral factor y = mi/ /mi. To 
be specific, for 6p = 7r/2 one finds an average transverse polarization of —0.3%, —3.7% 
and —0.02% for the three cases (/i — e), (r — > /i) and (r — > e), respectively. 

4 VF-boson propagator effects 

In order to incorporate IV-boson propagator effects one has to rewrite the charge re- 
tention form of the four-Fermi interaction (11) in terms of the charged currents of the 
Standard Model including the IV-boson propagator. This is easily done using the Fierz 
transformation property written down in Eq. (Bl) in Appendix B. 

For the present purposes it is only necessary to take into account ly-boson propagator 
effects in the Born term contributions. The momentum transfer is now q = — Pv^ = 
Pe + Pue- The q^q'^ piece of the W^-boson propagator contributes only at 0{mlm'j^/m^r) 
in the spectrum and rate functions and can therefore be dropped. In fact, using the Fierz 
identity (B2), it is not difficult to compute the contribution of the q'^q'^ piece exactly. 
The VT-boson propagator effect on the spectrum can thus be taken into account by the 
replacement 



\ mf^ x{2 - X 



m^J 3-2x 
^1 "^M 



m^y 1 — 2x 



1.3 , (24) 

2.4 . (25) 



where terms of 0{ml/m^r) have been neglected. Numerically the propagator corrections 
are quite small since m'i/m'^^ = 1.73 ■ 10~^ and ml/rn^^ = 4.88 • 10~^. 
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As is evident from Eq. (24) the H^-boson propagator affects tlie two pairs of spectrum 
functions differently. This means that one cannot absorb the IV-boson propagator effect 
entirely into a redefinition of Fermi's coupling constant Gp-, as advocated in Refs. [16, 17], 
when one considers polarization effects. 

In order to determine the propagator corrections to the rate functions one has to do 
the integrations 



/„ 



cixxGfril + ^^^^) = ±1(1 + -^), (26) 



dxxG^,r{l - -^^) = T^(l + ^^), (27) 
rriyy 1 — 2x 6 5 m^r 

where terms of 0{y) have been dropped in the Born term factors and in the integration 
measure. Again it is evident from Eq. (26) and Eq. (27) that the VT-boson propagator 
affects the two pairs of rate functions differently. If the VT-boson propagator effect is 
absorbed into a redefinition of Gi? using a measurement of the total unpolarized rate then 
this must be compensated for by multiplying the polarization dependent pieces G2,4 by 

2 

(1 — j:^) when calculating the average of the longitudinal polarization of the electron 

W 

(Pg) or the forward-backward asymmetry {Afs). 



5 0(a) corrections to spin dependent rate functions 

Many of the technical ingredients that go into the calculation of the 0{a) corrections can 
be found in a detailed account of the 0{as) corrections to the decay of a polarized top 
into a bottom quark and a W gauge boson t ^ b + presented in Ref. [2] (for the loop 
contribution see also [18]). Compared to Ref. [2] one needs to include the polarization 
dependence of the final massive fermion in this application. For the 0{a) tree-graph and 
one-loop contributions this is easily done. 

Let us briefly discuss the tree graph contribution. The 0(e) tree graph amplitude 
("internal bremsstrahlung" ) consists of the two contributions where the photon is either 
radiated off the electron or off the muon. One thus has 



where k and are the momentum and the polarization four-vector of the photon. Four- 
momentum conservation now reads = p^. + Q + k where Q is again the momentum 
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transferred to the neutrino pair. When squaring the tree graph amphtude one only sums 
over the polarization states of the photon since the muon and the electron are taken as 
polarized. Omitting again the antisymmetric contribution one obtains in the Feynman 
gauge 



1 (k-p^ + m^ Pe-pA^a-^ 



ik-Pi,)\ (k-p^) (k-pe)^ 



(k'^p^ + k^p'^ - k-p^g'^f^^ 



+ J^2(PePi + PeP", - Pe-P.g''^) ' J^^^P^e + P^Pe " P^'Pe^?"^) 



In the last line of Eq. (29) we have isolated the infrared singular piece of the charge- 
side tensor which is given by the usual soft photon factor multiplying the Born term 
contribution. In the phase space integration over the photon momentum the infrared 
singular piece is regularized by introducing a (small) photon mass which distorts the 
phase space away from the singular point. The infrared singular piece is cancelled by the 
corresponding singular piece in the one-loop contributions again calculated in Feynman 
gauge. In the one-loop amplitude the photon mass is included in the pole denominators 
only^. 

The remaining part of the charge-side tensor in Eq. (29) is infrared finite and can 
easily be integrated without a regulator photon mass. 

Just as in the Born term case treated in Sec. 3 the 0{a) charge-side tensor C^J is 
contracted with the neutral-side tensor A^°^ Eq. (16) where now Q = p^ — Pe — k. In 
the tree graph integration one first integrates over the photon's energy and then over the 
cosine of the angle between the photon and the electron taking appropriate care of the 
infrared singular piece. Any azimuthally dependent terms drop out after the azimuthal 
integration. One then finally adds in the one-loop contributions. The final result is 

®We specify our infrared regularization procedure since, historically, there has been a certain amount 
of controversy concerning the use of a photon mass regulator (see e.g. [19]). 
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TT TT izpx I ^ ' 



+ ((5 + 12a; - Slx^ + 28x^) - 3(19 - 40a; + ITx^)^^ - 3(19 - 4a;)y^ + Sy^jAa 

- 4/3x((ll - lOx + 5x2) _ 2(^1 + 53,)^2 ^ ^^^4\ 1 ^^^^^ 



^ ^^Bo™ ^ ^ -— ^(3/?V(l - 4x + 9y2) Ai - 8(l - X + y')' \ 
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TT ^ 7rl2/?2x2 

- (x(l + 33x2-28x^) - 3x(33-12x-17x2)y2 + 3(32 - 35x - Ax^)y^ - 5x/)A2 

- 2/3x((3 - lOx - 13x2 + 8x^) + (59 - llx^)y'^ - (35 - 10x)y^ + hy^) I, (31) 



^ «^Bo™ ^ ^ '^-^A - 3/5V(3 - 4x + 3y2)Ai - 8y2(l - X + y^)^ 
TT vr 12p^x^ [ ^ 

- (x(5 + 12x - 51x2 + 28x^) - (8 - 129x + GOx^ + 25x^)y2 

- 3(40 - 49x + 8x2)?/^ - (24 - 2?>x)y^ - 8y^)\2 



\3 



+ 2/3x((5 + lOx - 11x2 + 8x3) _ (^35 ^ i^^^^y^ + (59 _ ^Qx)?/^ + 3?/^) |, 
^ «^Bo™ ^ ^ ^^^1 - 3/?V((l - 4x)x + (8 + 3x)2/2) Ai 

TT TT 12p'^X'^ [ ^ ^ 

+ 8x(l - y2)(i _ X + 1/2)3^3 + 4/?x(x(l - 5x - 5x2 ^ 3^3^) 
- (4 - 59x + 14x2 + 5x3)?/2 - (104 - 59x + 5x2)y^ - (4 - x)y^ 
x2(l + 33x2 _ 283,3^) ^ (^4 ^ _ 1533,2 ^ iQ4^3 ^ 25x^)y2 



3(28 - 72x + 59x2 _ g^i^yi _ _ 24^ + 23x^)y^ + 4(1 + 2x)y^ A2 \ (33) 



(32) 
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^ »^Born s + «|/ _ (1 _ 2x + 3y') X, + ^i3x - 2x' - (4 - 3x)y')\2]. (34) 

We have used the abbreviations 

Ai = Iny, A2 = In , A3 = In ) ^ (35) 

and 



^ (l.>((l +,«)|) - i^)(A3 - A.) + (h{l - a.. 6)1) - l^)xl 



(36) 



The radiative corrections to the unpolarized spectrum function Gi including the full 
mass dependence have first been calculated in [20]. After correcting for an error in the 
calculation of [20], a correct result was published in [21]. Our results for G^"'' agree with 
those given in [21] and in [7]. Our results for Gg"'' agree with those of [7]. The results on 
G3 4 5 are new. 

Note that that the four spectrum functions ^^"234 are logarithmically mass divergent. 

We have checked that the expansion of the spectrum functions ^1^2 3 4 and G^^^ /y in terms 
of powers of y contains only even powers of y, in agreement with the general reasoning 
given in [22] . The leading term in the y-expansion of the spectrum functions Gi'2 3 4 can 
be reconstructed from the rrif, —>■ results to be presented in Sec. 6. We do not write 
down any of the higher order coefficients in the y-expansion, as was done in the Born 
term case, since such an expansion is not particularly illuminating. 

The explicit evaluation of the five spectrum functions is numerically quite stable except 
for the region very close to threshold, i.e. for /3-values smaller than f3 ^ 0.05, or, when 
expressed in terms of the scaled electron energy x, for x-values below (9.7 x 10~^; 1.2 x 
10~^;5.8 X 10~^). The origin of the instability are the inverse powers of /5 appearing in 
Eqs. (30-34). If one wants to explore the region below j3 ~ 0.05 one can make use of a 
small /3-expansion of the expressions Eqs. (30-36) which is not difficult to arrive at using 
an algebraic program such as e.g. MATHEMATICA. We do not write down explicit forms 
for the small /5-expansion because the expressions are not particularly illuminating. Let 
it be said that in the small /3-expansion G^i'} 5 and ^2^3 are even and odd functions of (3, 
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respectively. This shows that G2 and G3 are proportional to f3 just as in the Born term 
case. We mention that approximate formulas for the threshold region have been written 
down for Gi in [23]. 

In Fig. 2a (// —>■ e), Fig. 3a (r — * /i) and Fig. 4a (r e) we show plots of the x- 
dependence of the four spectrum functions jSxGi {i = 1, 2, 3, 4), with and without radiative 
corrections. The radiative corrections show a marked ^/-dependence. They are smallest 
for r ^ /i, become larger for /i — e, and are largest for r — » e. To a large part this can 
be traced to the (In?/) -dependent terms in the spectrum functions as will be discussed in 
more detail in Sec. 6. On an absolute scale the radiative corrections are generally quite 
small except for the hard end of the spectrum where they in fact (logarithmically) diverge 
^. On a relative scale the radiative corrections are quite large for (yU e) and for (r ^ e), 
and smaller for (r — /i) at the soft end of the spectrum, where the spectrum functions 
are small. This will show up in the radiative corrections to the longitudinal polarization 
of the daughter lepton P/, and the forward-backward asymmetry A pB which are large in 
the threshold region for fi e and r — > e and small for r ^ fi. As discussed in Sec. 3 
the contribution of (3xG^ is barely discernible for the case r — s> (Fig. 3a). At the scale 
of the figure the difference between the Born and NLO curves is not visible. 

It is interesting to note that the radiative corrections go through zeros close to x = 0.68 
and X = 0.82 for Gi^a and for G2,4, respectively, for all three cases discussed in this paper. 
The positions of the respective zeros are practically mass independent. Differences in 
the position of the zero show up only in the third digit. In fact, when discussing the 
mil — > case in Sec. 6 we have checked that the positions of the zeroes remain practically 
unchanged even when letting y 0. The radiative corrections are negative and positive 
below and above the zero for /5a;Gi^4, resp., and positive and negative below and above the 
zero for j3xG2^3- Qualitatively the alternating sign pattern over the range of the spectrum 
can be understood from the dominance of the (ln?/)-terms and from the fact that the 
(In I/) -dependent terms have to cancel out when one integrates over the spectrum. There 
is a tendency of the radiative corrections to cancel in the sums (Gi + G3) and (G2 + G4) 
and to add up in the differences (Gi — G3) and (G2 — Gi), i-e. the radiative corrections 
add destructively and constructively in the final electron's density matrix elements p++ 
and p , resp., to be discussed further in Sec. 6. 

We next turn to the radiative corrections of the longitudinal polarization of the 
electron and the forward-backward asymmetry Aps calculated according to Eqs. (9) and 
(20). We begin by discussing the limiting value of the longitudinal polarization at the 
soft end of the spectrum including the radiative corrections. At NLO one has 

^Thc divergent terms at the endpoint of the electron spectrum can be resummed into an exponential 
function [24, 25]. 



16 



lim p' = -Ipco^eJi-^ iji-yni + y') 1 

»-29 ' 3 \ ir72i/2+;(6-10t/-278i/'-10i/3+5i;'')-108;i^i/MriyJ 

For the cases fj, ^ e and r — e the term (IS?/^) in the second hne of Eq. (37) can be 
neglected and thus the hmiting value of the longitudinal polarization of the electron is 
given by Pj: ~ —-^P cos 9p which is smaller than the Born term value by a factor of 5. 
This is evident in Figs. 2b and 4b where the radiative corrections to the longitudinal 
polarization of the electron are shown. For the t —>■ fi case the term [ISy"^) dominates 
over (5/4)(a;/7r) and the limiting value of the longitudinal polarization of the /i is Pg ~ 
— ^PcosOp (1 — ^Ygp-), i.e. the correction to the Born term value is only PS 3.6%. This can 
be seen in Fig. 3b. The NLO limiting value for the forward-backward asymmetry is zero 
since, as discussed before, a small /5-expansion shows that G'^^ /G^^^ ~ {x — XminY^^ just 
as in the Born term case. This can again be seen in Figs. 2c - 4c. Finally, the behaviour 
of the transverse polarization P^ at x^m = 2?/ is quite similar to that of the longitudinal 
polarization. In fact, one just has to replace cos6p sinOp and 4 — > 5 in the numerator 
of the first line of Eq.(37) in order to obtain the limiting value P^. This leads to 

limP^^^-lpsm9pll--(^y^ + -] ]. (38) 

Just as for Pj: the 0{a) corrections are small for t —* fi and sizeable for fi ^ e and r — > e 
at threshold, as can be seen in Figs. 2b - 4b. Next we discuss the limiting behaviour of 
the spectrum functions at the hard end of the spectrum, where x x^ax = 1 + 1/^- As 
remarked on before the radiative correction contributions can be seen to logarithmically 
diverge in this limit (see footnote 7). Introducing x'^^^ = 1 + — e the limiting behaviour 
of the four spectrum functions is given by 



lim Gi = — lim G2 



a , 



lim G3 = lim 6*4 = (1 - vT + -2(1 - y ) x (39) 



vr 



X (1+y') T + ln 



Iny + (1 - y^) 1 + In 
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The limiting values of the four spectrum functions are, up to signs, all identical. This 
implies that the radiative corrections to the longitudinal polarization of the daughter 
lepton does not change the Born term value P// = — 1 at the hard end of the spectrum 
irrespective of the values of Pcos6p. Similarly the radiative corrections to the forward- 
backward asymmetry do not change the Born term value of Aj^b = — at the hard end 

of the spectrum. Finally, the NLO transverse polarization is zero at x = Xmax since ^5"^ 
is finite at Xmax- In fact, one has G^^^ —{a/7r)y{l — y^) \ny at Xmax- 

The fact that the relative corrections to the spectrum functions in the threshold region 
are relatively large for (/i e) and (r —>■ e) shows up in Figs. 2b and 4b, and in Figs. 2c 
and 4c, where the radiative corrections to the longitudinal polarization of the final state 
electrons and the forward-backward asymmetry are visibly large in the threshold region^°. 
Note, though, that we have enhanced the threshold region in our presentation of P// and 
ApB by choosing a logarithmic scale for x in Figs. 2b, c, 3b, c and 4b, c. For the forward- 
backward asymmetry the radiative corrections remain large over a larger part of the 
spectrum. The radiative corrections to P^ and ApB for the (r fi) decays are shown 
in Figs. 3b and 3c. As expected from the smallness of the radiative corrections to the 
spectrum functions shown in Fig. 3a the radiative corrections to both the longitudinal 
polarization of the final state muon and the forward-backward asymmetry are small. The 
same statement holds true for the radiative corrections to the transverse polarization. 
Since the spectrum functions are small at the soft end of the spectrum, the radiatively 
corrected average longitudinal polarization (P/,) of the daughter lepton /' is expected to 
remain very close to —1 in all three cases. Using the integrated rate functions presented 
at the end of this section we find the NLO result (P/,) = -0.999, -0.986 and -0.999 for 
{fj, — > e), (r — /i) and (r e), resp., for cos9p = with very little dependence on 
cos^p. The corresponding figures for (App) are —0.166, —0.170, —0.166. Finally, for the 
transverse polarization one finds (P^f ) = —0.0031, —0.037, —0.00018 for the three cases 
for cos^p = 0. 

Next we integrate the five spectrum functions f3xG\"^ over the electron spectrum ac- 
cording to (22). One obtains 



G^r^ = -1:^(1 - y'){75 - 956y2 + 75y^) - y\36 + y') y 

TT [ 4o 

vr^ 1 
- —(1 - 32?/3 + IQy^ - 32y^ + - -(60 + 270y^ - 4y^ + 17y^)y^ Iny 
4 6 

"'^'^In our numerical results for the longitudinal and transverse polarization of the daughter lepton and 

the forward-backward asymmetry we have not expanded out the inverse of the denominator function in 

terms of powers of a. 
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^(l-y^)(17-6V+17y^)ln(l-2/^) 



2(1 - y)\l + % + IQy' + %3 + y^) ln(l - y)\ny 



2(1 + y)^(l - % + lOy^ - %3 + y^) ln(l + y) In?/ 



+ (3 + 32y3 + 48y^ + 32^^ + 3y«)Li2(-2/) 



+ (3 - 32y3 + 48y^ - 32y5 + 3y«)Li2(y) 



^(«) ^^h^_y2^^^^y2^ ^3^4 _ 3^6)^1^(1 _ ^) ^ 2 ln(l + y)) \ny 

TT 6 



1 



432 



(1 - y)^(617 - 842?/ + 192%2 - 1592^^ - 3415y^ + 54y^ - 567y^ 



1 



- -(1 - y)y{l2 - 18y - 238^^ + 41y=^ - 1003y^ - 45?/^ - 92/^) Iny 
3o 



36 



(1 - 2/)2(13 + 261/ + 87y2 - 3641/3 + 535y^ - 102i/^ - Sly^) ln(l - y) 



+ Uj\U + 32y2 _ 3i/4) In^ y - 1(1 - y2)(25 + 13i/2 + 67i/^ + ISi/*^) ln(l + y) 
- Sy^i^- \v?y+ (2U^{-y) + U^^y) - Iny - 6Li3(-i/) - hh{y) - lc(3) 
+ 1(7 + 241/2 + 48i/^ - 8i/6 + 9/) ('Li2(-i/) + ^) |, 



a ( 1 



TT [3 

1 

432 



(1 - i/2)(3 - 131/2 _ y4 _ ^6)(i^(i _ ^) ^ 21n(l + y)) Ini/ 
(1 - i/)2(567 - 54y + 3415i/2 + 1592i/3 - 1929i/^ + 842i/^ - 6171/*^ 



(1 - i/)y(36 + 98y + 590i/2 - 265^^ + 27y^ - 99i/^ - 87y^) Ini/ 

36 
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- y)2(51 + 102?/ - 535?/=^ + 36V - 87i/^ - 2Qy^ - I2,y^) ln(l - y) 

36 

^7/2(8 + 66y2 - 2V - Z/') In' y - + 67y2 + iSy^ + 2hy^) ln(l + y) 

6 9 



7r 



7 



+ 8y^( - In^y + I 2Li2(-y) + Li2(y) - y 1 Iny - 6Li3(-y) - UM - -((3) 



+ (9 _ 82/2 + 482/^ + 2%*^ + 72/«) Li2(-2/) 



7r 



12 



(42) 



[432 



(l-y^)(581 + 61401/2 + 581?/^ 



- — (7 - 6y + 28y2 + 222y^ - 164y^ + 390y^ - 140y^ + 66y^ - 5y^) 
36 

- 4(1 - ?/')(13 - 176y2 + I3y4) ln(l - + 1^2(2 + 45^2 ^ 33^4 ^ ^6) 1^2 ^ 
36 3 

2 2 
+ -y{l - y)\5 -Ay + 5y^) ln(l -y)\ny- -y{l + y)\h + Ay + 5y^) ln(l + y) \ny 



-y2(212 + 930y2 + 388y^ - \ny 

-(1 - lOy - 562/2 - 102y^ - IQAy^ - 102i/^ - 56?/^ - lOy^ + y^)U2{-y) 
3 

^(1 + lOy - 56?/2 + I02t/^ - 164?/^ + 102?/^ - 56/ + lOy^ + y^)U2{y) |, (43) 



7r io 



^2/2(9 + - /) ln2 y + ^(1 - 32/2 + 32/ - 3/ + /) 
3 6 

^(1 - y)\l + Ay + 2/2) ln(l - 2/) In 2/ - ^(1 + 2/)'(l - 42/ + 2/') ln(l + y)\ny 

4(15 + 2492/2 + 1412/^ - 292/^) In 2/ + ^(1 - y^){ll + 382/2 + ll2/^) ln(l - y^) 
io 9 
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6G< 



/i - 

6Gr 



6G^ 



T 



6g: 



6G> 



SG,. 



6G< 



r — >• e 
5G> 



SGi 



+2.80% 
+8.22% 
+2.56% 
+7.82% 
+0.76% 



-2.69% 
-3.70% 
-2.68% 
-3.70% 
-3.68% 



-0.42% 
-0.68% 
-0.53% 
-0.79% 
-2.89% 



+0.76% 
+2.51% 
+0.66% 
+2.53% 
+0.09% 



-1.14% 
-1.59% 
-1.14% 
-1.59% 
-1.62% 



-0.40% 
-0.62% 
-0.45% 
-0.67% 
-1.49% 



- 5.20% 
-14.70% 

- 4.95% 
-14.28% 

- 1.72% 



-4.44% 
-6.10% 
-4.43% 
-6.10% 
-6.03% 



-0.42% 
-0.68% 
-0.54% 
-0.80% 
-4.53% 



Table 2: 

Numerical values of partially integrated {SGf , ^Gf) and total rate functions {SGi) divided 
by their respective Born term values. The symbols "<" and ">" stand for integrations from 
threshold to the zero point of the respective 0{a) contributions, and from the zero point to the 
endpoint of the spectrum. 



-(3 - 9y' - 32y' - 9y' + 3y%i2{~y) 



- o (3 - 9y' + 32y' - 9y' + 3y^)U2{y) 



(44) 



Note that G'^^ and G^^^ contain trilog functions, and associated with them, Euler's Zeta 
function C(3), whereas Gi^\ G^^^ and G^^^ contain only dilog functions. In agreement 
with the Lee-Nauenberg theorem [3] the rate functions do not contain any logarithmic 
mass singularities. Our result for G*^"^ agrees with the result in [26] where a different 
route of phase space integrations was taken to arrive at the total rate. Our result for 62"^ 



agrees with the result in [7]. The results on G^^\ Gi"'* and Gf' are new. As an additional 
check we have checked that all five rate functions G\"^ vanish for y ^ 1. 

In order to get a quantitative feeling about the size of the radiative corrections to 
the respective spectrum functions we have listed in Tab. 2 the percentage changes 6Gi 
induced by the radiative corrections, where 



Born 



Born. 



SG,, 



Born 



G 



(a) 



Born ' 



(45) 



Gi Gi 

The relative radiative corrections 6Gi can all be seen to be close to the naive expectation 
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of 0{a) where the relative radiative corrections to are largest. However, as was 



emphasized earlier on, all radiative corrections discussed in this paper go through zeros. 
Integrating over the whole spectrum therefore does not give an adequate representation 
of the size of the radiative corrections to the spectrum since there are sizable cancellation 
effects. This cancellation would become less effective if moments of the spectrum functions 
were taken. The moments could be chosen such that they either emphasize the threshold 
or the endpoint region. Alternatively, one can consider partially integrated rates where 
the integrations either run from threshold to the point where the radiative corrections go 
to zero or from the zero point to the endpoint. The two partially integrated rate functions 
will be denoted by Gf (lower part) and by Gf (upper part). The two respective partially 
integrated relative rate functions SGf and 6Gf are also listed in Tab. 2. The relative 
radiative corrections for the partially integrated rate functions can be seen to be much 
larger than for the fully integrated rate functions and can amount up to 0(10%) which 
is much larger than the naive 0{a) expectation. Note, though, that, in contrast to the 
rate functions, neither the moments of the spectral functions nor the partially integrated 
rates are free of logarithmic mass singularities. 

Since the complete rate expressions are rather unwieldy it is useful to consider the 
small ^/-expansions of the rate expressions. One has 




(46) 





(47) 





(48) 





(49) 
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a 



TT 



35 - Stt^ + 151ny 1 



(27 + 7r2 + 251n?/ + 



1 fi 

+ 12\n'y)y' + ^nY + 0{y') 



(50) 



It is well-known that the 0{a) small-y corrections to the reduced rate Gi"^ start only 
at 0{y^) (see e.g. [22, 27]). In contradistinction and contrary to the Born term case the 
mass corrections to the spin dependent rate functions G^2'\ G^^^ and G^^^ all start at 
0{y). 

In order to obtain a quantitative feeling about the importance of mass effects in the 
0{a) radiative contributions we have listed in Tab. 1 the percentage changes in the rate 
functions Gi when going from m;/ — > to m;/ 7^ using the m;/ —>■ results listed in Sec. 6. 
The mass effects are not small, in particular for the case (r — > /i). The percentage changes 
for the radiative contributions are larger than those in the Born term case which is partly 
due to the difference in the power pattern of the final state lepton mass corrections in the 
two cases. The quality of the m;/ approximation for the radiative corrections can be 
assessed by referring again to Tab. 1. The final state lepton mass effects tend to reduce 
the overall size of the radiative corrections. The reduction is largest (0(10%)) for the 
case (r — > /i) and smallest (O(10~^)) for the case (r e). They are largest for the rate 
functions G3 and G4 and smallest for Gi. That the mass effects are smallest for Gi is due 
to the fact that the mass corrections to Gi set in only at 0{y'^) (see Eq. (47)). Whether one 
is willing to tolerate the error incurred in using the simpler m// — » radiative correction 
formulas depends of course on the accuracy required for the application at hand. 

We now turn to the discussion of the NLO average longitudinal and transverse polar- 
ization of the electron {Pi) and (P^), and the NLO average forward-backward asymmetry 
(Afb)- They take remarkably simple forms in the y ^ limit. Including the Born term 
contribution and expanding the inverse denominator in powers of a/vr one has at NLO 

\nn{Pl) = -{l-^). (51) 

Note that (Pg) does not depend on PcosOp in this approximation. 
For the forward-backward asymmetry one obtains 

,• / . \ 1 r^/. " evr^ - 49, 

hmiApB) = --P(l ). (52) 

y^o b TT 9 

Finally, for the transverse polarization one obtains in the same approximation {cos 6 = 0) 
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, 1 , 2 , a 65 + 60 In , , 

lim(P,^) = - Py{l + ^^). (53) 

The 0(a) corrections to (Pe), {Afb) and to {Pf!~) are thus quite smalL The actual 
numerical values for (Pg), (Afb) and for {P^) listed earlier in this section lie very close 
to the above estimates in all three cases. 



6 The rUe — > limit and the anomalous helicity flip 
contribution 

The purpose of this section is two-fold. First we discuss the rrie limit of the mg 7^ 
0{a) results given in Sec. 5. This allows us to make contact with the mg results 
derived previously [6]. Second we discuss in some detail the origin of the anomalous 
helicity flip contribution resulting from collinear photon emission of the electron. Our 
results are presented in terms of the two diagonal components of the density matrix of 
the final state electron, which, in the limit mg 0, are nothing but the helicity no-flip 
and helicity flip contributions of the final state electron. The naive prediction of massless 
QED is that the helicity flip contributions vanish in all orders of perturbation theory. 
However, as first pointed out by Lee and Nauenberg [3], there will be a nonzero helicity 
flip contribution from collinear photon emission which survives the mg limit. This 
will be demonstrated in our 0{a) mg expressions. Our mg — > result for the helicity 
flip contribution is found to be in agreement with expectations derived from the universal 
equivalent particle approach of Falk and Sehgal [5]. For a discussion of the quality of the 
mi' approximation for the three cases (/x — e), (r — > /x) and (r e) we refer to the 
discussion at the end of Sec. 5. 

Since we want to discuss the helicity no-flip and flip contributions separately it is 
convenient to choose a slightly different representation for the differential rate Eq. (6). 
We write the differential rate in terms of the no-flip and flip contributions 

l-cos^) + 3-T---(l + cos^), (54) 



dx dcos 9p dx dcos 9p dx dcos 9 p 

where the no-flip and flip contributions are given by 

= _;3.r„((G. T G3) + (G. ^ G,)Pcos9p). (55) 

The terminology helicity no-flip {nf) and helicity flip [hf) is really only appropriate in the 
mg — i> limit where the electron emerging from the left-chiral weak interaction current is 
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purely left-handed. After photon emission the electron can then remain left-handed (nf) 
or can become right-handed (hf). For mg 7^ the respective nf and hf contributions are 
nothing but the (unnormalized) diagonal elements of the density matrix of the electron, 
i.e. T"-^ ~ p„„ and T^^ ~ p++. 

As concerns the helicity flip contribution one notes that there are no Born term helicity 
flip contributions in the limit rrig = since a massles electron emerging from the weak 
{V — A) vertex is left-handed. This is explicitly seen by inserting the Born term rate 
functions Eq. (19) in Eq. (55). Naively, one would expect no helicity flip contributions 
also at 0{a) because, in massless QED with rUe = 0, photon emission from the electron 
is helicity conserving. However, taking the limit rrie — >■ in Eqs. (30) - (33) one finds 
helicity flip contributions which survive the mg — limit. In fact, one finds 



a 



To({l-x)\5-2x) - {l-x)\l + 2x) PcosOp), (56) 

which agrees with the result presented in Ref. [6]. Because of the naive expectation 
that the helicity flip contribution vanishes in massless QED the presence of a helicity 
flip contribution is sometimes referred to as the anomalous helicity flip contribution. 
Moreover, the authors of Ref. [28] were able to show that the well-known axial anomaly 
can be traced to the existence of an anomalous helicity flip contribution to the absorptive 
part of the VVA triangle diagram in massless QED. This gives further justification for 
the use of the terminology "anomalous helicity flip contribution". In order to set the 
anomalous contributions apart we have highlighted them in Eq. (56) by enclosing them 
in square brackets. 

The authors of Ref. [6] had already expressed surprise at the simplicity of the struc- 
ture of the 0{a) helicity flip contributions without, however, attempting to identify the 
source of this structural simplicity. The simplicity of the helicity flip contribution be- 
comes manifest in the equivalent particle description of /i-decay where, in the peaking 
approximation, /^-decay is described by the two-stage process fi~ — ^ e~ followed by the 
branching process e~ — > e~ -|- 7 characterized by universal splitting functions Dnf/hf{z) 
[5]. In the splitting process z is the fractional energy of the emitted photon [5]. The 
off-shell electron in the propagator is replaced by an equivalent on-shell electron in the 
intermediate state. Since the helicity flip contribution arises entirely from the coUinear 
configuration it can be calculated in its entirety using the equivalent particle description. 

The helicity flip splitting function is given by Dhf{z) = az/ {27r), where z = ko/E' = 
[E' — E) / E' = 1 — x/x', and where ko is the energy of the emitted photon [5]. E' and E 
denote the energies of the initial and final electron in the splitting process. The helicity flip 
splitting function has to be folded with the appropriate mg = Born term contribution. 
The lower limit of the folding integration is determined by the soft photon point where 
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E' = E. The upper limit is determined by the maximal energy of the initial electron 
E' = One obtains 



dxdcosOp 2tt Jx x' dx'dcosOp 

= —To dx'{x' - x) ((3 - 2x') + (1 - 2x')P cos Op 
2% Jx ^ 



a 

12^' 



ToUl-xf{5 - 2x) - (1 - xf {l + 2x)P cos 9p 



(57) 



which exactly reproduces the result of Eq. (56). Note that the flip spectrum function does 
not contain a logarithmic mass factor. Integrating over the spectrum one obtains 



a. 



dcosOp IT 
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P cos 9p . 



{51 



which can be checked to agree with the results in Sec. 5 setting y = 0. 

The helicity no-flip contribution is again obtained by taking the rrie — >■ limit in 
Eqs. (30) - (33) but now for the differences of the respective spectrum functions as specified 
in Eq. (55). Including the 0{y^) Born term contributions one obtains 



dV'f 



dx dcosOp 



To x2(3-2x) + 



a 



(1 -x)2(5-2x) 



4x(ll - lOx + 5x^ 



- 6(6 - x)x ln(l - x) - 2(5 + 12a; - \5x^ + 4x^) In f ^ 



6(3 — AlX^x^ In 



\ — X 



12(3-2x)xV 



+ \x^{\-2x) 



a 



(l-x)2(l + 2x) 



2(3 - 10a; - I'ix^ + 



-(4 - 12x + 18x^ - 13x^) ln(l - x) + 2(1 + 21x^ - 4x^) In - 
X \x . 



6(1 -4a;)a;2 In 



1 — a; 



12(l-2x)a;V^ Pcos^p 



(59) 
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where 



a = 




1 — X 



X 




(60) 



The no-flip contribution agrees with the result presented in Ref. [6]. We have highlighted 
the anamolous contributions in Eq. (59) by enclosing them in square brackets. When 
calculating the total spectrum, i.e. when summing the flip and no-flip contributions (56) 
and (59), the anomalous contributions cancel. 

The 0{a) no-flip contribution is much larger than the flip contribution. This is illus- 
trated in Fig. 5 for the electron spectrum in the /i — *• e decay where the flip contribution 
has been multiplied by a factor of 20 in order to make it visible at all. The vanishing 
of the flip contribution for cos 9 = +1 at the hard end of the spectrum can be under- 
stood from angular momentum conservation. In order to be able to discuss the residual 
mass dependence the no-flip contribution has been split into its constant part and its 
logarithmic (In?/) part which come in with opposite signs over most of the spectrum, i.e. 
they partially cancel in the spectrum. Considering the numerical values of the ratios 
In / In ?/(^_^e) = 0.53 and Inyi^r^e)) / ^'^y{^l^e) = 1-53 it is clear from Fig. 5 that the 
cancellation between the constant and the logarithmic part is strongest for r — and 
weakest for r — > e. This observation provides a qualitative explanation of the hierarchy 
of the size of radiative corrections to the three decay cases as described in Sec. 5, namely 
the radiative corrections are largest for r ^ e and smallest for r ^ fi. 

We now turn to the discussion of the longitudinal polarization of the daughter lepton 
in the mg — > limit. It takes a rather simple form when one expands out the inverse 
denominator in terms of powers of {a/ir) keeping only the 0{a) contribution. In terms 
of the flip and no-flip contributions the longitudinal polarization of the electron reads 



When expanding the inverse denominator in terms of powers of (a/vr) one sees that the 
"normal" contributions in the numerator and denominator cancel exactly at 0{a) and 
one just remains with the anomalous contributions (No such cancellations occur for the 
forward-backward asymmetry ApB- We therefore refrain from presenting a closed formula 
for AfB in this approximation). One obtains 



pi 

e 



(61) 




(62) 
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It is clear that Eq. (62) does not apply very close to threshold. Due to the factor {l—xY/x^ 
in (62) the radiative corrections to the longitudinal polarization of the daughter lepton 
are largest close to threshold as is evidenced in the plots Fig. 2b and Fig. 4b. In the case 
of (r — > /x), mass effects prevent the radiative corrections to become large in the threshold 
region. 

Integrating the no-fiip contribution over the spectrum one obtains 



dcosOp \2 TT 



25 1 2\ 

H IT ] + 

16 4 / 




1 

24 



617 - 8Att^ 
432~ 



[Pcos^P . (63) 



which agrees with the results in Sec. 5. In particular one reproduces the simple expression 
Eq.(51) for the average of the longitudinal polarization. We have again highlighted the 
anomalous contributions by enclosing them in square brackets. The anomalous contri- 
butions can be seen to cancel in the spectrum and rate functions when adding up the 
spectrum no-fiip contribution (59) and flip contribution (56) , and the respective rate con- 
tributions (63) and (58). Numerically the 0{a) no-fiip contribution dominates over the 
0{a) (anomalous) flip contribution. For the 0{a) contributions to the rate functions, 
which do not depend on y in the y ^ approximation, one obtains (—0.168, —0.121, 
-0.107) for the rate ratio dr'^-^("V^^r"^^"^ for cos^p = 1,0, -1. 



7 Summary and conclusions 

We have computed the 0{a) corrections to the leptonic decays of the /x- and r-leptons 
including polarization effects and the full mass dependence of the respective final-state 
leptons. The radiative corrections to the spectrum functions are sizeable for (/i e)- 
decays, large for (r e)-decays, and smaller for (r yu)-decays. In large part this 
pattern is due to the (In i/)-dependent contributions to the spectrum. The polarization 
of the final-state lepton deviates substantially from the naive m// = values P// = — 1 
and P;f = towards the soft end of the spectrum. The radiative corrections to the 
longitudinal and transverse polarization of the daughter lepton in the threshold region 
are substantial for (/i e)- and (r — > e)-decays and small for (r /i)-decays. Similar 
statements hold for the forward-backward asymmetry ApB- 

For the rate functions we have compared our 0{a) rriii ^ results with m// results 
derived previously in [6]. In particular in the (r fi) case the errors incurred in using 
the 0{a) mil — > results are large (of 0(10%) in the 0{a) rate functions). A mass effect 
is already showing up in the experimental values for the branching ratios of the two decay 
modes t~ ^ fi~ + + Vj. and r~ — > e~ -|- z/g + Vt- They are BR{t~ ^ fi~ + + v^-) = 
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(17.37 ± 0.06)% and BR{t~ ^ e" + i^^ + i^r) = (17.84 ± 0.06)% [16]. The two branching 
ratios are compatible with the mass dependence of the Born term rates. In order to be 
sensitive to the mass dependence of the radiative corrections the error on the branching 
ratios would have to be improved by at least a factor of ten. 

Whether one is willing to tolerate the error brought about by using the simpler m^/ 
radiative correction formulas depends of course on the accuracy required for the applica- 
tion at hand. We nevertheless strongly recommend use of the complete results in numerical 
investigations instead of using the mii — > approximation. The analytical m;/ 7^ formu- 
las written down in this paper are of sufficient simplicity to allow for easy incorporation 
into numerical programs. 

From what was being said in Sec. 3 it is clear that the results of this paper can im- 
mediately be applied to the case of semileptonic quark decays, where e.g. in the case of 
the semileptonic 6 — > c + /" + z// decays the final state c-quark mass can certainly not be 
neglected. 
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A Appendix A 

When integrating the spectrum functions it is convenient to transform to the integration 
variable ^, where x = 7/(^2 + l)/{, such that P = {1 - ^^)/(l + ^^). 
Trilog functions are generated from the integrals 

^\ny+\\n^y-Usiy) + aS), (Al) 
6 3 

-\n{y)Lhiy')-Lk{y)+U,{y'), (A2) 

^ Hy)U2{y') + 2 Us{y) - U^iy') - C(3), (A3) 

^ Iny - i Hy)U2{y') - 2 U,{y) + U,{y') + C(3), (A4) 
6 z 
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y 1 

1 e 



y 1 



!-!/{, 



where 

LW.) := - /' dy, Li3(.) /■ L!eM (A5) 

Jo y JO y 

Use has been made of the relation 

Li„(x) + Li„(-x) = ^^Li„(x^). (A6) 
Euler's Zeta function is defined by 

oo 

C(5) = E^"^ C(3) = 1.202057... (AT) 

k=l 

B Appendix B 

The Fierz identity 

- 1^)Vp [7m(1 - = - [7m(1 - 15)U (Bl) 

is well known. Not so well known is the Fierz identity (see e.g. [10]) 

[(1 ± ^,)U [(1 T = ^[7^(1 T 75)],/? [7m(1 ± 75)]... (B2) 

Latter identity allows one to transform the q'^q'^ piece of the VT-boson propagator dis- 
cussed in Sec. 4 back into the standard form N'^^Cap used in the remaining part of the 
paper. 
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Figures 2a, b, c: 

The case ;U — > e. Scaled energy dependence of 

a) spectrum functions P xGi, i = 1, . . . ,5 

b) the longitudinal polarization of the electron for cos 9p = —1, 0, +1, and 

c) the forward-backward asymmetry Aps for cos 9 = 0. 

All curves with and without radiative corrections. Polarisation P is set to P = 1. 
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Figure 5: 

Scaled energy dependence of flip (hf) and no-flip (n/) spectrum functions in the limit y 
0. Latter contribution is separated into its non-logarithmic part (n/;non-lny) and its mass 
dependent logarithmic part {nf; Iny) which is plotted for the case (// — > e). 
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